Study of °He+ C Elastic Scattering Using a Microscopic Optical Potential 
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The ^He+^'^C elastic scattering data at beam energies of 3, 38.3 and 41.6 MeV/nucleon are studied 
utilizing the microscopic optical potentials obtained by a double-folding procedure and also by using 
those inherent in the high-energy approximation. The calculated optical potentials are based on 
the neutron and proton density distributions of colliding nuclei established in an appropriate model 
for ®He and obtained from the electron scattering form factors for ^^C. The depths of the real and 
imaginary parts of the microscopic optical potentials are considered as fitting parameters. At low 
energy the volume optical potentials reproduce sufficiently well the experimental data. At higher 
energies, generally, additional surface terms having form of a derivative of the imaginary part of the 
microscopic optical potential are needed. The problem of ambiguity of adjusted optical potentials is 
resolved requiring the respective volume integrals to obey the determined dependence on the collision 
energy. Estimations of the Pauli blocking effects on the optical potentials and cross sections are also 
given and discussed. Conclusions on the role of the aforesaid effects and on the mechanism of the 
considered processes are made. 

PACS numbers: 24.10. Ht. 25.60.-t, 21.30.-x. 21.10.Gv 



I. INTRODUCTION 

Experimental and theoretical studies of exotic light nu- 
clei with a localized nuclear core and a dilute few-neutron 
halo or skin have been an important and advanced area 
in the nuclear physics in the past decades. The availabil- 
ity of radioactive ion beams facilities made it possible to 
carry out many experiments and to get more information 
regarding the structure of these nuclei and the respective 
reaction mechanisms (for more information see, e.g., the 
recent review of the problem in Ref. [l|). In this sense, 
^He is a typical nucleus having the weak binding energy 
and extended neutron halo in its periphery. The latter is 
the reason why in collisions with the proton and nuclear 
targets the projectile nucleus ^He is breaking up with a 
comparably large probability that causes the flux loss in 
the elastic channel. Therefore, the study of elastic scat- 
tering of ^He on protons or light targets is a powerful 
tool to get information on peculiarities of the mechanism 
of such processes. 

The data on cross sections of processes with light ex- 
otic nuclei have been analyzed using various phenomeno- 
logical and microscopic methods. Among the latter 
we should mention the microscopic analysis using the 
coordinate-space g-matrix folding method (e.g., Ref. 0] 
and references therein), as well as works where the real 
part of the optical potential (ReOP) is microscopically 
calculated (e.g., Ref. [31]) using the folding approach (e.g., 
Refs. 0-0|)- Usually the imaginary part of the OP's 
(ImOP) and the spin-orbit terms have been determined 
phenomenologically that has led to the usage of a num- 
ber of fitting parameters. In our previous works (8|, |9| in- 
stead of using a phenomenological imaginary part of OP 
we have performed calculations of ^He-|-p 8] and ^He+p 
elastic differential cross sections by means of the mi- 
croscopic OP with the imaginary part taken from the 



OP derived in [T^, U^ in the frameworks of the high- 
energy approximation (HEA) |12h14| that is known as 
the Glauber theory. In the case of ^He-|-p elastic scat- 
tering, it has been shown in our previous study [S] that 
the depth of the imaginary part of the respective micro- 
scopically calculated OP has to be appreciably changed 
to get an agreement with the existing experimental data. 
The present study of ^He-|-^^C scattering could give a 
novel information on mechanism of the process due to 
the more complicated dependence of the microscopic OP 
not only on the density of the projectile ^He but also on 
the density of the target nucleus. 

In the last years a number of works has been devoted 
also to the elastic scattering of ^He on ^^C nucleus and, 
particularly, to the study of the mechanism of this process 
including the role of breakup channels. In the present pa- 
per we perform an analysis of the ^He-I-^^C elastic scat- 
tering data at three beam energies E = 3 [l^, 38.3 |la | 
and 41.6 MeV/nucleon [131 using the microscopically cal- 
culated OP. The data have been already considered indi- 
vidually in the frameworks of other theoretical models. In 
Ref. [15] the differential cross sections at total energy 18 
MeV were analyzed by means of Woods-Saxon (WS) OP 
with radius parameter of the imaginary part Rj — 5.38 
fm being about twice larger than that of the real part 
{Rj^ — 2.4 — 2.7 fm). The so large absorbtion radius 
in elastic channel may thought to be caused by breakup 
channels which take place at the far periphery of OP. 
Recently, the same ^He-|-"'^^C elastic scattering data were 
fitted in Ref. [18] by OP having a squared WS real part 
and a standard WS shape for its imaginary part. Fur- 
thermore, a larger radius Ri ~ 6.17 fm was obtained for 
the latter. The different forms of the real part of OP and 
also the values of Rj in Rcfs. [15] and 18] just reflect the 
known problem of the ambiguity of parameters of phe- 
nomenological OP's when fitting them to the restricted 



amount of experimental data. In principle, this problem 
does not arise in microscopic OP's whose basic param- 
eters have been already established by fitting to other 
data. In this line, a part of the problem was overcome in 
Ref. ^i, where the elastic scattering data at E ~ 38.3 
MeV/nucleon were analyzed using the semi-microscopic 
OP. Its real part was defined in a double-folding model 
including the direct and exchange convolution integrals, 
while the imaginary part was taken phenomenologically 
in the WS form. Adjusting the latter to the data at rela- 
tively small angles, the reduced radius and diffuseness pa- 
rameters were obtained in the range of r/ = 1.471 — 1.569 
fm and aj = 0.358—0.524 fm. Then, to get a better agree- 
ment at larger angles, the dynamic polarization potential 
(DPP) in the form of a derivative of a WS function was 
added to the volume OP. It affects strongly the total OP 
in the surface region at radii around 4—5 fm. 

Recently, in Ref. [19| five Gaussian-like forms for 
the ^He matter density distributions were tested using 
the real part of OP without the exchange term [J| to- 
gether with the 5-parameter phenomenological imagi- 
nary part (volume and surface) based on the WS form. 
This model was adopted to study the data at 38.3 and 
41.6 MeV/nucleon. Also, the authors generated another 
(microscopic) OP by involving in folding calculations 
the complex Jeukenne-Lejeunne-Mahaux (JLM) effective 
nucleon-nucleon (NN) potential [2^. Doing so, the op- 
timized values of the fitted parameters were fixed and 
then, to improve slightly the fits to data the repulsive 
real DPP term with two free parameters was introduced, 
as well. In this way, the ambiguity problem is retained 
just in such semi-microscopic models of OP. 

In Ref. |21| , an attention was paid to the dynamic po- 
larization term in OP. A Monte Carlo method was de- 
veloped to calculate both the ^He variational wave func- 
tion and the Glauber amplitude of the microscopic scat- 
tering of the nucleons of ^He by the nucleus ^^C as a 
whole. To this end the phenomenological p+^'^C OP was 
utilized and, as a result, the full ^He-|-^^C OP was re- 
stored from the calculated Glauber eikonal phase. Then, 
the difference between this OP and the single-folding OP 
estimated without accounting for the Glauber multiple- 
scattering terms, was defined as DPP responsible for the 
breakup channels. One can see from Fig. 5 of Ref. [21[ 
that at 40 MeV/nucleon in the range of 4— 5 fm the imag- 
inary single-folding potential (negative) W is about 11 — 3 
MeV and the DPP makes it deeper by about 7—3 MeV, 
so the full eikonal W has to be about 18—6 MeV. The cer- 
tain way to investigate simultaneously the effects of elas- 
tic and breakup processes can be performed in the frame- 
work of coupled-channel (CC) methods. In Refs. |13, ^^ 
these methods have been elaborated to estimate the ef- 
fects of the elastic and inelastic breakups of the projec- 
tile nucleus ^He— > a+n+n on the elastic scattering of 
^He-I-^^C. Two models for the ^He structure were ex- 
plored in Ref. 17]. One of them uses the modified wave 
function from [23] of the 3-body a+n+n system con- 
sidered as a "Borromean" one. The second model con- 



structs the a- "di-neutron" potential and the correspond- 
ing 2-body wave function. So far the input potentials 
Ua+„, Ua+2n and C/q+i2c, C/„+i2c, [/2„+i2c havc been 
taken from the respective fits and estimations. These 
both models are completely parameter-free and explain 
the elastic scattering data at 41.6 MeV/nucleon fairly 
well. The analysis of the breakup effects of ^He on the 
elastic scattering is done in Ref. [22,] using the so-called 
continuum-discretized CC method. The reaction system 
is described as a four-body model of a+n+n+^'^C The 
wave function of the bound and continuum states of the 
3-body system a+n+n is presented applying the speci- 
fied Gaussian expansion method, and the ground state 
wave function of ^^C is calculated by the microscopic 3a 
cluster model. The resulting microscopic OP was calcu- 
lated as double-folding integral with the calculated den- 
sities of ^He and ^^C, and then multiplied by the complex 
factor [Nb, + iNj) with coefficients optimized by a fit to 
the elastic scattering data. The exchange terms in the 
folding OP were neglected. The results for the ^Ke+^^C 
scattering showed that the optimum value of Nj is 0.5 
at 3 MeV/nucleon and 0.3 at 38.3 MeV/nucleon, respec- 
tively, while Nn = 1. Also, it was shown that the ef- 
fects of the coupling of channels are more important at 
comparatively large angles. The coupling smooths the 
diffraction structure of the differential cross sections at 
E — 3 McV/nucleon and shifts down the corresponding 
curve at i<^ = 38.3 McV/nucleon calculated without cou- 
pling. 

It can be mentioned that calculations by coupled reac- 
tion channel models with accounting for the cluster and 
continuum states are encouraged to study their sensitiv- 
ity to the input information on the reaction mechanism. 
On the other hand, the breakup reactions reveal them- 
selves through the dynamic polarization terms in the full 
OP for elastic scattering. The explicit information on 
these channels can be obtained from the unambiguous 
OP obtained from the respective analysis of the elastic 
scattering experimental data. In our paper we try to re- 
alize the following idea. We start analyzing the elastic 
scattering data by the microscopic optical potential ob- 
tained in Ref. [l3]. Its real part includes the direct and 
exchange terms that are the same used in Ref. [laj. The 
imaginary part of OP is based on the Glauber theory of 
high-energy scattering of complex systems and is an inte- 
gral which folds the nucleon-nucleon scattering amplitude 
/nn with the density distribution functions of the bare 
nucleons of colliding nuclei. Therefore, first, this OP con- 
sists of parameter-free real and imaginary parts defined 
by the respective terms of Jnn- However, this poten- 
tial reveals only the single-particle physical nature of the 
colliding nuclei because it depends on the single-particle 
nuclear densities. Then, we assume that additional terms 
to our basic OP (which have to be added to explain the 
experimental data) may be considered as a consequence 
of the presence of more complicated channels. In the 
case of the loosely bound ^He projectile these terms are 
thought to arise due to the breakup channels. Thus the 



main effort should be directed to minimize the ambigui- 
ties in the fitted OP's by studying differential elastic cross 
sections at different energies and to involve external phys- 
ical conditions in order to make as narrow as possible the 
corridor of the deviations of selected theoretical curves. 
The paper is organized as follows. The theoretical 
scheme to calculate the real and imaginary parts of the 
OP is given in Sec. |lTl The results of the calculations of 
OP's and the ^He-|-^^C elastic scattering differential cross 
sections and their discussion are presented in Sec. IIIII 
The summary of the work and conclusions are given in 
SecHVl 



II. THE MICROSCOPIC OPTICAL POTENTIAL 

Here we give the main expressions for the real and 
imaginary parts of the nucleus-nucleus OP 



and with terms of the energy and density dependencies: 



U{r) ^V'^''{r)+iW{r). 



(1) 



The real part V^^ consists of the direct and exchange 
double-folding (DF) integrals that include an effective 
NN potential and density distribution functions of col- 
liding nuclei: 



V°^{r) = V"ir)+V^^ir). 



(2) 



The formalism of the DF potentials is described in de- 
tails, e.g., in Refs. i^iSj]. In general, in Eq. ^ V^ and 
V^'^ are composed from the isoscalar (IS) and isovector 
(IV) contributions, but in the considered case the isovec- 
tor part is omitted because Z = N in the target nucleus 
^^C and, thus, one can write: 

y^(r) = j d\pd\tpj,{vp)pt{vt)v^j,{s), (3) 

V^^{r) = d^rpd^rtppirp,rp + s)ptirt,rt-s) 

iK(r) • s 



xi;f^(s)exp 



M 



(4) 



where s = r 4- r^ — r^ is the vector between two nucleons, 
one of which belongs to the projectile and another one to 
the target nucleus. In Eq. ([3]) Ppi^p) and pt(rt) are the 
densities of the projectile and the target, respectively. In 
Eq. dH) pp{vp,Yp + s) and pt{rt,rt — s) are the density 
matrices for the projectile and the target that are usu- 
ally taken in an approximate form [24] . In the modern 
calculations of the DF potentials the effective interaction 
v^^ (of CDM3Y6-type) based on the Paris NN forces 
and having the form 

v^j^iE,p,s)^g{E)F{p)vis) (5) 

is usually applied with the distance dependence given by 






(6) 



g{E) = 1 - 0.003E;, F{p) = C 



1 + ae 



-lip 



IP 



(7) 



The energy dependent factor in Eq. ([7]) is taken to be 
a linear function of the bombarding energy per nucleon, 
while p in F{p) is the sum of the projectile and target 
densities, p = pp + pt- The parameters Ni, fXi [Eq. (j6|)], 
C, a, (3, 7 [Eq. ([7])], and all details of the mathematical 
treatments and calculations are given in Refs. (g. [25| . 



In Eq. Q 



EX 



is the exchange part of the effective 



NN interaction. It is important to note that the energy 
dependence of V^^ arises primarily from the contribu- 
tion of the exponent in the integrand of Eq. (j4]). Indeed, 
there the local nucleus-nucleus momentum 
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l/^^(r)-K(r)] 



1/2 



(8) 



with Ap, At, m being the projectile and target atomic 
numbers and the nucleon mass, and M = ApAt/{Ap + 
At). As can be seen, K{r) depends on the folding poten- 
tial V^^{r) that has to be calculated itself and, there- 
fore, we have to deal with a typical non-linear problem. 
Usually, two different kinds of effective NN potentials are 
employed in calculations, namely the Paris CDM3Y6 and 
the Reid DDM3Y1 NN interactions, which are defined 
by two different sets of the aforementioned parameters. 
The direct parts of these potentials have different signs 
and, for example, in the case of CDM3Y6 forces the V^^ 
is negative while V^ is positive. So, if in the calculations 
one takes only the direct part of V^^ with the Paris M3Y 
NN forces, then the corresponding real part of such OP 
is positive one. Therefore, one should proceed carefully 
when neglecting the exchange part of OP. 

Concerning the imaginary part of our OP, we take it 
in two forms. In the first case the imaginary part has the 
same form as the real one but with different strength. 
At the same time we test another shape of the imagi- 
nary part that corresponds to the full microscopic OP 
derived in Refs. (TSUll within the HEA [Ulil. In the 
momentum representation this OP has the form 

F 1 

^''^''pp{q)pt{q)fN{q)d\. (9) 

Here a^ and aisr are the averaged over isospins of nuclei 
the NN total scattering cross section and the ratio of 
real to imaginary parts of the forward NN amplitude, 
both being parameterized, e.g., in Refs. [2y, [23|. The 
NN form factor is taken as /Ar(g) = exp(— g^/3/2) with 
the slope parameter /3 — 0.219 fm^ [2^. It is easy to 
verify that the eikonal phase for this potential 



i=l 



1 f 



z^)dz 



(10) 



is reduced to the HE A microscopic phase ^n for nucleus- 
nucleus scattering. To this end, let us substitute ^ in 
(fTO| and simplify the latter in the cylindrical coordinate 
system where d^q = q±dq±d(l)dq\\ and qr — q±b cos ip + 
Q\\z, q^ = ql +gj|, q± = qcos{d/2), q^ = gsin(i9/2). 
Here 1} is the angle of scattering, b is the impact parame- 
ter of the projectile trajectory of motion directed straight 
ahead along the z axis. In HEA one neglects the longi- 
tudinal part of the transfer momentum q^^ ^ q± and the 
small angle terms q\\z 2± (7||2i? ~ kRi}^ <C 1 to give q ~ 
q± at 'd <^ ^Jl/kR. Therefore, the smooth functions be- 
come p{q) — >■ p[qi_), and fN{q) — >■ fNiq±), and then one 
can perform integrations J_ dq\\ exp(— ig||z) = 2TrS{z) 

and /p dcj) eiq){—iq±b cos <p) — 2'KJo{q±b). Thus, as a re- 
sult we obtain the standard form of the HEA phase in 
the so-called optical limit of the Glauber theory: 



47r 



X / Jo{qb)pp{q)ptiq)fNiq)qdq. (11) 

Here one sets q± — q — 2ksm{t')/2) that is valid at small 
angles of scattering ■& <^ ^Jl/kR. 

As a rule this phase is employed to estimate the HEA 
elastic scattering amplitude 



f{d) ^ ik 



Mqb) 



1 



J^N{h) 



bdb. (12) 



In applications of f{'d) from Eq. (fT2|) the main limita- 
tions E :s> \U\ and 1} <C ^1/kR are connected with the 
basic assumption that the integration in Eqs. (|10l) - p^ 
is performed along the z-axis with a straight-line classi- 
cal trajectory of motion. To correct partly this approxi- 
mation at lower energies and larger angles one can take 
into account the trajectory distortion. For the Coulomb 
distortion the respective prescription was used, e.g., in 
Ref. [231 , where the impact parameter b in ([12]) was re- 
placed by be = a-|-Va^ + b"^ with a = ZiZ2e'^ /2Ec.m. be- 
ing the half-distance of closest approach in the Coulomb 
field of a point charge. The approximated method for the 
inclusion of a distortion in presence of a nuclear potential 
was formulated, e.g., in Ref. [30|- However, all these at- 
tempts can not fully account for distortion effects of the 
classical trajectory caused by the total complex OP. In- 
stead, the conventional way to resolve this problem is to 
compute numerically the Schrodinger equation with the 
initial OP given by Eq. ([9]). Moreover, when using this 
original OP in the wave equation, it becomes not neces- 
sary to neglect the longitudinal terms in the momentum 
transfer. In this way, one expects that the initial OP 
([9]) can be used not only for the scattering at high ener- 
gies but also for comparably lower energies and for wider 
range of scattering angles. 

Hereafter we shall use only the imaginary part of 
the full OP © transformed (by using the equality 



/ di^q exp (— iqr) = ATrjo{qr)) to the form 



X / Joiqr)pp{q)ptiq)fN{q)q dq. (13) 
Jo 

In the further calculations the microscopic volume op- 
tical potential has the following form: 

U{r) = NnV^^ir) + iNiW{r), (14) 

where W{r) is taken to be equal either to V'^^{r) or to 
W^{r). The parameters Nji and Nj entering Eq. P^ 
renormalize the strength of OP and are fitted by compar- 
ison with the experimental cross sections. In the present 
work, attempting to simulate the surface effects caused 
by the polarization potential |3ll - [33| , we add to the vol- 
ume potential [Eq. p^ ] the respective surface terms. 
Usually, they are taken as a derivative of the imaginary 
part of OP, as follows: 



W^ir) = -iN\ 



sf 



= -i.N' 



s/. 



dWir) 
dr 
dW{r) 



^iN'/r 



dr 
,dW{r) 



dr 



= -iN' 



,fdW(r-S) 



dr 



(15) 
(16) 
(17) 
(18) 



where Nj is also a fitting parameter, S gives the shift of 
the potential ([T8|) and in our case is fixed to be (5 = 1 fm. 
Concluding this section, we would like to emphasize 
that our basic OP contains the same real part as that 
one in Ref. [1^, but instead of the phenomenological 
ImOP, we utilize a microscopically calculated imaginary 
potential. In addition, in contrast to Ref. [l9| where 
only the direct part of the ReOP and a part of the 
ImOP were calculated microscopically, in the present 
work we include also microscopic exchange part of the 
OP. Concerning the comparison with the experimental 
data, we consider not only selected data (as e.g., at 
E = 38.3 MeV/nucleon in iS] , and at £; = 38.3 and 41.6 
MeV/nucleon in [19]) but add also the data at fairly low 
energy E — 3 MeV/nucleon and analyze simultaneously 
the three sets of data. This allows us to determine the 
ambiguities when adjusting the values of the OP param- 
eters [Eqs. p^ - (IT5|) ] because we include an additional 
condition in the fitting procedure. Namely, we consider 
also the energy behavior of the volume integrals of ReOP 
and ImOP as an additional physical constraint. In this 
way, the information on the dynamical polarization po- 
tentials obtained from such more precise analysis can be 
considered as more reliable. 



III. RESULTS AND DISCUSSION 

In this section we present the results of the calcula- 
tions of the microscopic OP's and the respective elastic 
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FIG. 1: Upper part: total, point-proton and point-neutron 
microscopic LSSM densities of ^He (from Ref Q) (a) and the 
density of ^^C [H.lli] (b). Lower part: microscopic OP's l/-"^ 
(c) and W" (d) for the "^He+^^C elastic scattering at S = 3, 
38.3 and 41.6 MeV/nucleon (iV^ = Ni ^ 1 and iV;-'' = 0). 



scattering differential cross sections at energies E < 100 
MeV/nucleon obtained following the theoretical scheme 
in Secini In contrast to the cases of ^'^He-l-p elastic scat- 
tering where only the density of the projectile ^'^He takes 
part in the calculations, in our case both densities, of the 
projectile ^He and the target ^^C, have to be included in 
the calculations of the OP (see Eqs. ^ and (g])). The 
results of our work Q on the ^He-|-p elastic scattering 
at i5 < 100 MeV/nucleon showed that the large-scale 
shell model (LSSM) density of ®He microscopically calcu- 
lated in Ref. ^] using a complete Ahuj shell-model space 
and the Woods-Saxon single-particle wave function ba- 
sis with realistic exponential asymptotic behavior is the 
most preferable one and it is used also in the present 
work. For ^^C we use the symmetrized Fermi- type den- 
sity with the radius and difFuseness parameters c = 3.593 
fm and a — 0.493 fm from Ref. [34]. They were obtained 
by defolding the ^^C charge density distribution deduced 
in Ref. [35] from analysis of the corresponding electron 
scattering form factors. 

In Fig. [1] are shown the densities of ^He and ^^C, as 
weh as the OP's V"^ calculated using Eqs. ^-^ and 
W" obtained within the HE A [Eq. ^] for the three 
cases of the incident energy that are considered {E — 3, 
38.3 and 41.6 MeV/nucleon). It can be seen that the 
increase of the energy leads to reduced depths and slopes 
of ReOP and ImOP. 

We calculated the ^He-I-^^C elastic scattering differen- 
tial cross sections using the program DWUCK4 [3^ and 
the microscopically calculated OP [Eq. P^ ]. As already 
mentioned in Sec. |TT1 in the calculations we add to these 



volume potentials the respective surface terms [Eqs. ([15])- 
^TE\i ]. The latter is done only for the ImOP, having in 
mind that due to the breakup channel effects there is 
a "loss of the flux" from the elastic channel. We note 
that for the ^He-t-^^C process there is not a spin-orbit 
contribution to the OP in contrast to the ^^^He-t-p cases. 

In our work we consider the set of the Ni coefficients 
{Nn, Nj and Nj) as parameters to be found out from 
the comparison with the empirical data. We should men- 
tion (as it had been emphasized in our previous works 
[3, Q for ^'^He+p scattering) that we do not aim to find 
a complete agreement with the data. The introduction of 
the A^ 's related to the depths of the different components 
of the OP's can be considered as a way to introduce a 
quantitative measure of the deviations of the predictions 
of our method from the reality (e.g., the differences of 
A^'s from unity for given energies). 

The starting energy of our calculations was E — 38.3 
MeV/nuclcon. At this energy HEA can be applied as a 
good approximation to calculate the ImOP. As a first ex- 
ample, we present in Fig.[2]the results of our calculations 
using: i) only volume OP P^ for the two types of ImOP 
{V^^ and W^) and ii) different forms of the surface con- 
tributions to the ImOP [Eqs. (fT5t - (fT8)) ]. The parameters 
A^'s are determined by a fitting procedure. The results 
of the calculations are close to each other and that is 
why all of them are presented inside a gray area. Two 
definitions of x^ are used: 
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(19) 



(20) 



In the first definition the x^ values were obtained consid- 
ering uniform 10% errors for all the analyzed data. This 
procedure is often used by other authors. Here (j*'^{di) 
and a'^^'P{'di) are the theoretical and experimental values 
of the differential cross sections {da/dfl) or their ratio to 
the Rutherford cross section. In the last case the xi i^ 
dimensionless and its values for the results in Fig.[2]range 
in the interval 0.191 < xi — 0.362, while the values of 
Nr are in the interval 0.893 < A^ < 1.268. The values 
of the predicted total reaction cross section a^ are also 
calculated. 

One can see from Fig. [5] that the inclusion of differ- 
ent forms of the surface potential leads to almost similar 
results for the cross section. This was also the case of 
^He-|-p processes studied in our previous work [9]. As 
is known, the problem of the ambiguity of the values of 
A's arises when the fitting procedure concerns a limited 
number of experimental data. 

The situation is ambiguous also for the case of the 
energy E = Al.6 MeV/nucleon. 

The case oi E = 3 MeV/nucleon is a particular one be- 
cause of this rather low energy. Nevertheless, we made an 



10 



D 
^ 1 



10" 



11 1 1 1 1 1 1 1 
; E=38.3AMeV 


; 


• t 






-. 


+ 
, ... 1 ... . 


■ 



10 15 

9c.m. [deg] 



20 25 



a/ 



FIG. 2: Cross sections of He+ C elastic scattering at E 
38.3 MeV/nucleon calculated by fitting the Nr, Ni, N^ 
parameters of the microscopic OP [Eqs. p4|l - (ll8t ] (gray area) 
The experimental data are taken from Ref. [l6l |. 



attempt to consider it using OP obtained in our method. 
The calculations showed that for this energy the fitting 
procedure led to Nj = for the case of the surface term 
given by Eq. ([T5|) . We note that in the case of _E = 3 
MeV/nucleon the ambiguity in the explanation of the 
data [1^ still remains. 

In what follows, we tried to choose the most physical 
values of iV's for the energies considered. As is known, 
the fitting procedure belongs to the class of the ill-posed 
problems (e.g., Ref. |33). To resolve this problem it is 
necessary to impose some physical constraints when fit- 
ting the parameters of a model. In our case it might be 
the data on the total cross sections but often the corre- 
sponding values are missing. Another physical criterion 
that has to be imposed is the obtained potentials to obey 
a determined behavior of the volume integrals |^] 



Jv 
Jw 



An 
ApAt 

An 
ApAt 



NnV^'^'irydr, (21) 

NiW{rydr, (22) 



as functions of the energy. Indeed, it was shown for 
nucleon and light-ions scattering on nuclei (see, e.g., 
Refs. |38l - l4q |) that the values of the volume integrals 
Jv decrease with the energy increase at £■ < 100 
MeV/nucleon, while Jw increases at low energies up to 
10-20 MeV/nucleon and then saturates. We would like 
to note that such conditions were also imposed in Ref. \9\ 
when the microscopic OP's were introduced to study the 
^He+p scattering and their depth parameters Nr and Nj 
were fitted. The values of Jy and Jw for the ^He-f-^^C 
scattering that fulfil the condition for their energy depen- 
dence are presented in Tables HI [ll] and IIIIl In the cases 
when we include surface terms to the ImOP we modify 
Jw accounting for them. 

In the next part of the work we select those sets of the 
parameters A^'s that lead to the already mentioned be- 
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FIG. 3; The energy dependence of the volume integrals Jv 
and Jw that are related to the selected OP's from a number 
of fitted microscopic OP's with and without surface terms. 
The values of N's, Jv, Jw, X^ ''■nd an corresponding to the 
symbols in panels (a), (b), and (c) are given in Tables HI HTl 
and mil respectively. 



havior of Jy and Jw as functions of the energy for three 
cases: 1) only the volume terms; 2) the volume terms 
plus the surface term given by Eq. (jl6|) : 3) the volume 
terms plus the surface term from Eq. (J17p . The behav- 
ior of the volume integrals as functions of the energy is 
presented in Fig. |3l 

Using the same values of the parameters A^'s already 
selected, we present in Figs. HIE] and H] the cross sections 
for the three energies and for the three cases mentioned 
above. One can see that the best agreement with the 
data for all the three energies can be obtained by the OP 
with the volume and surface term [Eq. (1171) ] whose vol- 
ume integrals follow (though approximately) the already 
mentioned energy dependence. 

In Figs. [Tlja) and [TJ^b) are given those microscopic 
ReOP and ImOP that lead to the best agreement with 
the experimental data for the ^He-|-"'^^C elastic scatter- 
ing cross sections shown in Figs. SJ [S] and [B] for ener- 
gies E = 5, 38.3 and 41.6 MeV/nucleon. For £: = 3 
MeV/nucleon there are only volume OP's, while for 
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FIG. 4: Differential cross section of elastic He+ C scatter- 
ing at _E = 3 (a), 38.3 (b) and 41.6 MeV/nucleon (c) calcu- 
lated using only volume OP [Eq. (fTi)) ]. Solid line: W = W'^ , 
dashed line: W = V^^ . The values of the fitted parameters 
Nu and Ni corresponding to the curves in the upper, middle 
and lower part are given in Table |I1 The experimental data 
are taken from Refs. I15l-[l7ll. 
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FIG. 5: The same as in Fig. |4] (without the case for E = 3 
MeV/nucleon) at E ^ 38.3 (a) and E = 41.6 MeV/nucleon 
(b) calculated using volume OP and surface contribution to 
the ImOP [Eq. (|16p ]. The values of the parameters A^'s are 
given in Table [III The experimental data are taken from 
Refs. Hill. 



TABLE II: The same as Table |I] but for the parameters Nu, 
Ni, and Nj of the total OP with the surface term from 
Eq. ^. 



E 



W Nr Ni N 



s/ 



Jv 



Jw 



CTR 



38.3 W" 1.000 0.023 0.082 278.740 109.172 17.399 1055.67 

38.3 V'^^ 0.924 0.082 0.101 257.556 106.420 5.006 1174.66 

41.6 W" 0.852 0.337 0.051 232.645 188.590 3.734 1070.77 

41.6 y^^ 0.800 0.500 0.014 218.446 147.876 3.781 1072.40 



TABLE I: The optimal values of the parameters Nr, Nj for 
the volume OP [Eq. (fTi)) ] for the elastic ^He+^^C cross sec- 
tions at energies E — 3, 38.3 and 41.6 MeV/nucleon when 
the imaginary potential W was selected in the forms W^ or 
yDF rpj^g values of the volume integrals Jv and Jw, X^ ^nd 
total reaction cross sections aR (in mb) are also given. 



TABLE III: The same as Table U but for the parameters Nr 

Ni, and 

Eq. {TT)). The values of Nj^ are in fm'^ 



Nj of the total OP with the surface term from 



E 



W Nr Ni 



Jv 



Jw 



X 



O-R 



3 
3 
38.3 



W" 0.826 0.154 297.109 

V°^ 0.793 0.345 285.239 

W" 1.268 0.511 353.442 

38.3 V°^ 1.123 0.472 313.025 

41.6 W" 0.897 0.689 244.933 

41.6 V°^ 0.814 0.584 222.269 



212.952 9.121 1427.33 
124.095 9.890 1428.52 
208.567 80.808 1028.77 
131.565 50.847 1033.79 
265.680 3.737 1067.32 
159.466 3.774 1067.55 



W Nr Ni N; 



sf 



Jv 



Jw 



o-fl 



3 W" 
3 V°^ 
38.3 W" 
38.3 V°^ 
41.6 W" 
41.6 V°^ 



0.790 0.074 0.002 
0.725 0.040 0.008 
0.932 0.028 0.019 
0.932 0.204 0.012 
0.797 0.255 0.011 
0.578 0.041 0.022 



284.160 137.506 
260.779 55.924 
259.786 110.017 
259.786 105.469 
217.627 152.281 
157.827 98.546 



8.912 1449.98 

9.418 1533.22 

5.059 1185.22 

8.425 1161.91 

3.711 1091.46 

2.398 1224.91 
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FIG. 6: The same as in Fig. 4 but for the volume OP and sur- 
face contribution to the ImOP [Eq. (|17|) ]. The corresponding 
A'^'s are given in Table UTIl The experimental data are taken 
fromRefs. [IMS. 



E ~ 38.3 MeV/nucleon there is a surface contribution 
to the ImOP using W = V^^ in Eq. ^ a.nd W = W" 
in Eq. (ITTI) : for E ~ 41.6 MeV/nucleon there are only 
volume OP's. 

In Figs.jSJb) and[5l^c) the real and imaginary parts of 
the OP for _E = 3 MeV/nucleon obtained in the present 
work are compared with the phenomenological OP's from 
Ref. [1^ (where WS forms have been used for ReOP and 
ImOP) and from Ref. fl8] (with OP having a squared 
WS real part and a standard WS shape for the ImOP). 
The results for the cross sections are shown in Fig. |8l[a) . 
One can see much better agreement for our cross sections 
obtained using microscopic OP's than those obtained in 
a phenomenological way in Refs. [l5l[l8|. 

We should note also that in Ref. ^ the ReOP (Vo) 
increases and the ImOP (Wq) decreases with the energy 
increase which is in contradiction with the generally ac- 
cepted results and with the behavior of the volume inte- 
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FIG. 7: Selected OP's that lead to best agreement with the 
data of the ®He-|-'^'^C cross sections shown in Figs. 21 O and 
|6l The values of the parameters Af's for the curves at _E = 3 
and 41.6 MeV/nucleon (with VK = W^ and W = V°^) are 
listed in Table [11 while at i? = 38.3 MeV/nucleon are given in 
TableHfor W = V°^ and Table [Till for W = W^ . 



grals as functions of the energy [38H40| . 

In Fig. [3] the deviations of the OP's from their volume 
parts are presented for the three energies considered. In 
our theoretical scheme they are related to the surface 
parts of the ImOP in the form of Nj rdW^''^^^ /dr and 

Nj r^dW'^^^^^ /dr. As is known, these contributions 
can be considered as the so-called dynamical polarization 
potential that owes its origin to effects of the breakup of 
a pair of neutrons from ^He. 

In what follows, we would like to discuss the deviations 
of the values of the coefficients N's (and, particularly, 
of TV/) from unity. As known, a folding model can be 
thought as meaningful only when the renormalization co- 
efficients of the folded potential are close to unity. In our 
case (see Table [J) Ni = 0.154 for .B = 3 MeV/nucleon 
and Ni = 0.689 for E = 41.6 MeV/nucleon. Here we 
would like to emphasize that the obtained values of Ni 
within the HEA for the small energy E ^ 2> MeV/nucleon 
reflect the effects of Pauli blocking (see, e.g. [41]) that in 
the case of nucleus-nucleus scattering reduce the depth 
of ImOP by a factor of 10 or more. This is related to the 
fact that in the HEA the microscopic optical potential is 
proportional to the free NN cross section (ctat), that in 
nuclear matter is reduced by the so-called in-medium fac- 
tor fm that accounts for, in particular, the Pauli blocking 
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FIG. 8: (a) Differential cross section of elastic ^He+^'^C scat- 
tering at _E = 3 MeV/nucleon. Solid and dotted lines show the 
results with microscopic ImOP W'^ and V^^ , respectively. 
The results with the phenomenological OP's from Refs. [l5|| 
and [18] are given by dashed and dash-dotted lines, corre- 
spondingly. The experimental data are taken from Ref. [15J : 
The ReOP and ImOP for iJ = 3 MeV/nucleon microscopi- 
cally obtained in the present work and those from Refs. [1^ 
and [ig] are given in panels (b) and (c), respectively. 
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There are many estimations of the factor fm, e.g. those 
based in the Briickner-Hartree-Fock (BHF) theory [43 - 
Hl]. In [43] the expression for fm is obtained by a least 
squares fit to the experimental total reaction cross section 
data over a wide incident energy range. The latter gives 
the separate forms for the Upp and crp„ cross sections 
accompanied by the factors 
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FIG. 9: The surface ImOP's (dynamical polarization OP's) 
used in the calculations of the cross sections of "He-I-^^C elas- 
tic scattering at i? = 3 (a), 38.3 (b), and 41.6 MeV/nucleon 
(c). Solid and dashed lines: using Eq. (|16|l with values of A'^'s 



from Tableland with W ^ W" and W = V"^^ , respec- 
tively. Dot-dashed and dashed two-dots lines: using Eq. (|17p 
with iV's from Table [Till and with W = W" and W = V°^ , 
respectively. 



fminn) 



- 7.772 £;0-06p'** 
1 + 18.01 pi-46 ■' 



(25) 



where E is the kinetic energy in laboratory system per 
nuclcon of the projectile nucleus, and p — pp + pt- (One 
can guess that in Eqs. (|24|) and ([25]) the numerical co- 
efficients in the second terms of fm have the respective 
dimensions that to measure Eiab in MeV and p in fni~~f). 
Recently, expressions for fm were presented in Ref. |44l | 
in the approximation of the isotropic free TV TV cross sec- 
tion. In this approach the Pauli projection operator in 
the G-matrix as a solution of the Bethe-Goldstone equa- 
tion is considered as a geometrical factor that restricts 
the conditions on the available angles of scattering of the 
NN pair to unoccupied final states. The result for the 
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NN correction factor is in the form [44j : 

1 



fn 



1 + 1.892 



f \pp-pt\ \ 

\ PPO J 



2.75 



f{E), 



where 



f{E > 46.27p2/3) = 1 - 



37.02 



f{E < A6.27p^'-^) 



■^/3\ - 



E 



231.38p2/3' 



(26) 



(27) 



(28) 



E is the laboratory energy per nucleon in MeV, po — 0.17 
fm^'' and p — [pp + pt)/ po. In practice, hard numerical 
problems arise when one intends to use these formulae in 
calculations of folding integrals for the microscopic po- 
tentials. Instead, it is easy to estimate the in-medium 
effects in the realistic case suggesting that the main con- 
tribution comes from the region of half-density radii of 
the colliding nuclei, where p = Pp + pt = Po, Pp ~ Pt = 
and p = 1. Then, the first term in Eq. ()26|) equals to 1, 
and thus one gets f{E) ~ 0.013 for E = 3 MeV/nucleon 
and f{E) ~ 0.18 for E = 41.6 McV/nucleon. The use of 
Eqs. dHJ and ([25]) (from (Hi) lead to the following esti- 
mations: fm{np) = 0.717 and fm{nn) = 0.68 for i? = 3 
MeV/nucleon and fm{np) = 0.75 and fm{nn) = 0.74 for 
E — 41.6 MeV/nucleon. One can see that the results for 
Ni = 0.154 (at ^ = 3 MeV/nucleon) and Ni = 0.689 
(at E = 41.6 MeV/nucleon) from our work mentioned 
above and obtained by a fitting procedure are in a cor- 
rect "direction" and they are between the estimations 
using Refs. H!] and [ii]. 

We note that these estimations are valid only for the 
volume OP's (even in the local density approximation). 
We emphasize that one has to account also for the com- 
petition with the channels at the nuclear periphery (the 
breakup) that, according to the coupled-channel calcu- 
lations, play an important role. Their contribution that 
had been initially obtained only accounting for the chan- 
nels of NN scattering inside the nuclear matter leads 
to changes of the ImOP in the elastic channel. An ex- 
ample for such a " renormalization" of ImOP for elas- 
tic d -\- A scattering due to the stripping channel was 
given in Ref. [45|. For more complex systems this is 
difficult to be done because only in the {d,p) reaction 
one can use the approximation of the "delta" -potential 
in the deuteron. Thus, we note that not only the Pauli 
blocking, but also other processes due to different mech- 
anisms (breakup and others) play role, as it is considered 
in various references mentioned in the Introduction of the 
present work. 



IV. SUMMARY AND CONCLUSIONS 

The results of the present work can be summarized as 
follows: 



(i) The microscopic optical potential and cross sections 
of ^He-f ^^C elastic scattering were calculated at the en- 
ergies of -E = 3, 38.3 and 41.6 MeV/nucleon. Compar- 
isons with the experimental data and results of other ap- 
proaches were presented. The direct and exchange parts 
of the real OP {V^^) were calculated microscopically us- 
ing the double-folding procedure and density-dependent 
M3Y (of CDM3Y6-type) effective interaction based on 
the Paris nucleon-nucleon potential. The imaginary parts 
of the OP were taken in the forms of V^^ or W^ , the 
latter being calculated using the high-energy approxima- 
tion. The microscopic densities of protons and neutrons 
in ^He calculated within the large-scale shell model were 
used. The nucleon density distribution functions of ^^C 
were taken as defolded charged densities obtained from 
the best fit to the experimental form factors from electron 
elastic scattering on ^^C. In this way, in contrast to the 
phenomenological and semi-microscopic models we deal 
with fully microscopic approach as a physical ground to 
account for the single-particle structure of the colliding 
nuclei. The elastic scattering differential cross sections 
were calculated using the program DWUCK4. 

(ii) While at low energies the volume OP's can repro- 
duce sufficiently well the experimental data, at higher 
energies additional surface terms in OP having a form 
of a derivative of the imaginary part of the OP became 
necessary and were used in the present work. 

(iii) The depths of the real and imaginary parts of the 
microscopic OP's are considered as fitting parameters. 
As is expected when one utilizes the fitting procedure 
in the case of a limited number of experimental data, 
the problem of the ambiguity of these parameters arises. 
To overcome (at least partly) this ambiguity, additional 
physical constraints should be imposed. Doing so, we 
require in our work the values of the depth's parame- 
ters A^'s to lead to volume integrals Jy and Jw with 
realistic energy dependence for energies < i? < 100 
MeV/nucleon. Namely, Jy's must decrease while Jh^'s 
increase to some constant values with the increase of the 
energy. 

(iv) The comparison of our results with those of some 
phenomenological approaches pointed out the advantages 
of using microscopic real and imaginary parts of the op- 
tical potential imposing realistic physical constraints on 
their depths as that one of the behavior of the volume 
integrals as functions of the energy. 

(v) As in works of other authors (e.g., Ref. [5|) we 
consider in more details the behavior of the OP in the 
nuclear periphery. This gives a possibility to make some 
conclusions about the contributions of the dynamical po- 
larization terms of the OP or, in other words, about the 
coupled-channel effects. 

(vi) It is shown that the deviations of the values of A^'s 
from unity (given in TableUfor the volume OP's) that are 
obtained by a fit to the experimental data for ^He-|-^^C 
elastic scattering are related to the Pauli blocking ef- 
fects. These values are smaller than unity and turn out 
to be between the approximate estimations performed 



using the results of approaches where Pauh blocking is 
taken into account (e.g. Refs. [ij, [ijl). It is also men- 
tioned that together with the important Pauli blocking 
effects, the role of other mechanisms, such as breakup 
processes, also have to be accounted for in the study of 
the reaction considered. 
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